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~"1 Abstract 

*V We demonstrate that the charge distributions in Hubbard-model representations of transition metal oxide heterojucntions can be 
described by a Thomas-Fermi theory in which the energy is approximated as the sum of the electrostatic energy and the uniform 
CZ5 three-dimensional Hubbard model energy per site at the local density equals to a constant. When charged atomic layers in the oxides 
are approximated as two-dimensional sheets with uniform charge density, the electrostatic energy is simply evaluated. We find that 



this Thomas-Fermi theory can reproduce results obtained from full Hartree-Fock theory for various different heterostructures. We 
also show explicitly how Thomas-Fermi theory can be used to estimate some key properties qualitatively. 
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, Recent advances in techniques for layer-by-layer growth 
of transition metal oxides promise new types of heterostruc- 
tures between strongly-correlated materials |H2l3j and 
has motivated many theoretical studies devoted to explor- 
ing their physical properties [4I5I6| using simple Hubbard 
models augmented by a long-range Coulombic interaction 
term. The charge distribution near a the transition metal 
oxide heterostructures (TMOHs) is usually studied using 
Hartree-Fock theory or dynamical mean-field theory with 
satisfactory results. We have found that the same results 
can be reproduced much more simply using a Thomas- 
Fermi theory (TFT) in which all energy contributions other 
than the electrostatic one are approximated by a local 
density approximation [516] . In this paper, we derive TFT 
specifically for the generalized Hubbard models of TMOHs 
from general consideration and demonstrate how it can be 
used to estimate some key properties of the TMOHs. 

The simplest toy model to investigate the TMOHs is: 
H = Hjj+Hc, where Htj is the single-band Hubbard model 
describing on-site correlations and He is the long-ranged 
Coulomb interaction. Following the spirit of TFT, we can 
write down TF equation as: 



ju(p[f]) + vnir) = const. 



(1) 



where p(p) and ujj(r) are the chemical potential of the 
Hubbard model for a uniform system with density p and 



the electrostatic potential at position r. Since the transi- 
tion metal oxides have a natural layered structures, we al- 
low the average density to vary from layer to layer near the 
heterojunction. It is usually a good approximation to ig- 
nore the periodic but non-uniform density variation within 
each layer in calculating the electrostatic potential. In this 
approxmation both u(f) and vh{t) depend only on the po- 
sition along the growth direction z and the TF equation 
reduces to: 



u(p(z)) + vh(z) = const. 



(2) 
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where p(z) is the average electron density at z. 

The electrostatic potential vr can be evaluated from 
self-consitently determined densities by solving the Poisson 
equation. In the continuum limit each layer can be con- 
sidered as a two-dimensional sheet with uniform effective 
charge density a. Consequently, the electrostatic potential 
at z is contributed by all charged layers at z 1 with effective 
charge density cr(z'): 

v h(z) = =2irU c J2<r(z')\z-z'\ (3) 

z' z' 

where U c = e 2 /ea is the strength of long-ranged Coulomb 
interaction, and a is the lattice constant. The next thing 
is to evaluate the chemical potential u(p(z)). Ideally 
u(p(z)) should be obtained from exact results for the 
homogeneous three-dimensional Hubbard model. Since 
this is not available some approximations are required. 
The simpliest choice is the Hartree-Fock approximation. 



Since Hartree-Fock theory can have separate solutions 
for antiferromagnetic (AFM) and ferromagnetic (FM) 
states(among others), we can solve TF equation for both 
cases. It has been shown that TFT accurately repro- 
duces the charge distribution along the growth direction 
obtained by Hartree-Fock theory for modulation-doped 
Mott-insulator-Mott-insulator and Mott-insulator-band- 
insulator heterostructures[5]. 

Recently the authors proposed that electronic inter- 
face reconstruction (EIR), a charge tranfer from the out- 
most surface layer to the interface, can occur in a polar 
(AMO3 )-nonpolar (A'M'O s ) peroskite Mott-insulator 
heterostructure[6]. The effective charge densities are a = 1 
for the AO layers, for A'O layers, p(z) for the MO2 
layers, and (p(z) — 1) for M'Oz, where p(z) is the average 
electron density on the <i-band of transitional metal ions 
M(M'). We consider a film with N p pairs of AO ~ M0 2 
layers and (N — N p ) pairs of A'O - M'0 2 layers. The 
electrostatic potential for this case is [6]: 



2nU r 



J2\z~z A \~y2 W) - 0(z' - N p )] \z' - z|(4) 
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where z A = 1.5, 2.5, ...,N P + 0.5, (z, z') = 1, 2, ...,N, and 
9(z' — N p ) is the Heaviside step function. The charge distri- 
bution calculated from TFT is given in Fig[TJ and is nearly 
identical to the results of Hartree-Fock theory presented in 
the Fig. 2(b) of Ref[B]. The success of TFT for different 
TMOHs is remarkable. The critical reason for this success 
is that the charge distribution in these TMOHs is a conse- 
quence of competition between the on-site correlation en- 
ergy (represented by p(p(z))) and the spatial long-ranged 
Coulomb field (represented by vh{z)), and these two dom- 
inant energy scales are correctly captured by TFT. This 
feature also promises that if the chemical potential p(p(z)) 
can be obtained more realistically by some sophisticated 
approaches for the homogeneous three-dimensional Hub- 
bard model, TFT can also predict the correct charge dis- 
tribution for most TMOHs without heavy computation. 

Another advantage of TFT emerges if we want to es- 
timate some key properties qualitatively. It is known, for 
example 0, that a minimum value of N p is required for EIR 
to occur. This minimum value can be estimated easily us- 
ing TFT. Since only the outer most surface and the inter- 
face layer have electron densites different from 1 when EIR 
first occurs, we consider: 

(J>Ui(ft(zi)) + v H {zi) = pu 2 (p( z i)) + vn(zi) (5) 



where Z\ = 1 is the position of outmost layer and Zi — 
N p + 1 is that of doped interfacial layer by our convention. 
U\ (U2) is the Hubbard U parameter for polar (nonpolar) 
material. We put p{z{) = 1 — x, pizi) = 1 + x, and p{z) = 1 
for all other values of z, then we can evaluate vh{z\) and 
vii(zi): 
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Fig. 1. Electron distribution calculated by TFT for AFM and FM 
states with Ui/t = U 2 /t = 20, U c = 0.8, N p = 5, and N = 10 for a 
polar (AM 03)-nonpolar (A 1 M' O3) heterostructure. 



vh(zi) 
2ttU c 

VH(zj) 
2TlU r 



N p 

— ^ \zi - z A \ - \z' - Zi\ - x\z { - z\\ 

Z A z'=2 

N p 

= \ z i - z a\ - e i z ' ~ Zi \ ~ (- 1 ~ x ^ zi ~ zi i 



(6) 
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It is straightforward to compute: ^2 ZA 



zi~z A \ = E Z4 
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, and 
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~ 2 ' t—iz'=2 \" *J-I — L-iz'=2 \" "*\ ~ 2 

z\\ = N p . So if we require x > 0, we find N p should 

be: N p > MLr2(1+ 2 t"tj7 l(1) - Thc insi g ht revealed by this 
expression has been discussed in Ref [6] . A similar analyis 
has been discussed earlier [5]. 

In summary, we have derived a TFT appropriate for Hub- 
bard model descriptions of TMOHs and applied it to cubic 
peroskite heterojunction systems. We have shown that this 
TFT can be used in an appropriate way to estimate some 
key properties which are difficult to treat analytically in 
Hartree-Fock theory and dynamical mean-field theory. The 
applicability of TFT indicates that the on-site correlation 
and long-ranged Coulomb interaction are the most domi- 
nant energy scales in TMOHs. For this reason it is possi- 
ble to create new types of two-dimensional systems which 
inherit strong-correlation features from the Mott insulator 
near interface between TMOs. 
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